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*~i I The Littlest Higgs with T-parity Model (LET) is one of the attractive candidates of 

physics beyond the Standard Model. In this study we focus on heavy gauge bosons 
and heavy leptons that LHT imposes, which are detectable at ILC. The mass and cou- 
plings of these particles possess important information of the model's mass generation 
mechanism and how they interact with other particles. Precision measurement will be 
essential among verifying the model. We will address the measurement accuracy of 
i-Q , Little Higgs particle masses, couplings and model parameters at ILC. 
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The Standard Model has made remarkable success in accurately describing particle physics. 
However the mass of the Higgs boson suffers from an instability under radiative corrections. 
This "hierarchy problem" leads us to expect new physics near the weak scale. In the Little 
Higgs scenario, the Higgs boson is regarded as a pseudo Nambu-Goldstone boson arising 
^Q ■ from the spontaneous breaking of a global symmetry. Due to a imposed symmetry, new 

[■■--. I particles such as heavy gauge bosons and top-partners are introduced. Contributions of 

^D ■ these particles cancel out the main quadratically divergent corrections of the Higgs boson 

^^ I mass at a one-loop level. Thus solving the hierarchy problem. 

f^ ■ In this study we discuss the Littlest Higgs model and furthermore implement a Z2 sym- 

Cn I metry called T-parity, which was imposed to satisfy constraints from electroweak precision 

measurements [iHS]. The lightest T-odd particle (heavy photon) is stable and provides a 
good candidate for dark matter. In order to test the Little Higgs with T-Parity (LHT) 
I^ . model, precise determinations of properties of little Higgs partners are essential, because 

V^ I these particles hold information of mass generation mechanism and thus are directly related 

C^ . to the cancellation of quadratically divergent corrections to the Higgs mass term. The mea- 

surements of heavy gauge boson masses and heavy leptons are particularly important. Heavy 
gauge bosons acquire mass through the breaking of the global symmetry and can solely be 
described with the parameter /, i.e. the energy scale of the global symmetry breaking. 
Heavy leptons also acquire mass terms through the breaking of the global symmetry and 
are described with / and the lepton yukawa coupling(K). All the masses and interactions in 
the gauge boson sector and lepton sector are described with these two parameters. We can 
extract the two mass originating parameters by measuring LHT masses. Furthermore LHT 
particles have characteristic interactions similar to the standard model, and they can also 
be explained with these two parameters. 

Here we will discuss the precision measurement of three heavy gauge bosons and two 
heavy leptons. Checking the consistency with extracted parameters from measured masses 
and couplings will be strong evidence that the discovered new particles are indeed LHT 
particles. 
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At the Large Hadron Collider (LHC), although it has potential of discovering LHT 
particles, because there are missing particles in the final state, the center mass energy is 
unknown and there are large background events, it is extremely difficult to reconstruct 
the masses of T-parity odd LHT particles. On the other hand, the ILC with it's clean 
environment and fixed center mass energy, will be an ideal environment to measure the 
properties of the LHT particles. 

In this study we will address ILC's measurement sensitivity to masses, parameters, couplings 
of heavy gauge bosons and heavy leptons based on realistic Monte Carlo simulation. 

2 Representative point and target mode 



In order to perform a numerical simulation at ILC, we need to choose a representative point 
in the parameter space of the LHT. Firstly, the model parameters should satisfy the current 
electroweak precision data. It also has to satisfy the cosmological observation of dark mat- 
ter relics. Thus, we consider not only the electroweak precision measurements but also the 
WMAP observation to choose a point in the parameter space. The current measurement of 
relic density is f2/)j\//i^= 1.06±0.008 from this measurement, we have selected a representa- 
tive point for the Higgs mass and vacuum expectation value (/), where they are 134 GeV 
and 580 GeV, respectively. 

As for the lepton Yukawa coupling(K), we have selected a representative point of k=0.5. 
If the yukawa coupling is too small heavy leptons become light and should be observed 
at previous experiments. If the yukawa coupling is too large, the four fermion interaction 
contribution increases and large discrepancy from the standard model emerges. With the 
parameter sets of {mn, /, K)=(134GeV, 580GeV, 0.5), the masses of the heavy gauge bosons 
and leptons are {Mah,Mwh,Mzh,M^hMuh) = (81.9 GeV, 368 GeV, 369 GeV, 410GeV, 
400GeV), where Ah^Zhi Wh, gh, i^h are the Little Higgs partners of a photon, Z boson, W 
boson, electron, neutrino respectively. Ah is a dark matter candidate. Since all the partners 
are lighter than 500 GeV, it is possible to generate them at ILC. 







v) v„v„ 



Figure 1: Feynman diagrams 



To measure heavy gauge boson properties we used three processes; e+e — > AhZh, ZhZh 
and WhWh- Ah Ah is also producible, but because Ah is a dark matter candidate the sig- 
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nal is undetectable. To measure heavy lepton properties we used two other processes; e"'"e~ 

— >■ grGh and vri^h- ^hZh was studied at yG = 500 GeV and the other processes at yG = 

1 TeV. The heavy gauge bosons Zh and W//100%branching ratios. Zh decays 100% to Ah 

and higgs, Wh decays 100% to Ah and W. Thus we have chosen the following final states 

for heavy gauge bosons. 

i) AhZh -^ AnAnh -^ AnAnbh 

ii) ZhZh -^ AnAnhh -^ AnAHqqqq 

iii) WhWh -> AhAhWW -^ AnAnqqqq 

At the representative point we have chosen, heavy Leptons are heavier than heavy gauge 

bosons, therefore they can decay emitting heavy gauge bosons. Branching ratios of ch and 

vh are shown at Table [TJ 



Table 1: Branching ratio of heavy leptons 



Particle 


Branch 


Branching ratio 


Particle 


Branch 


Branching ratio 




cAh 


30% 


I'M 


eWn 


65% 


en 


eZn 


25% 


vZu 


11% 




eWn 


45% 


vAh 


24% 



For the heavy leptons we have chosen the following final states, 
iv) eneu — > eeZuZn — > eeAnAnqqqq 
v) vhvh — > eeWnWH — > bcAhAhWW 

Among the three possible ways a heavy electron can decay, the i^Wh branch has difficulty 
accessing the mass of ch since en directly emits a missing neutrino and a darkmatter emit- 
ting Wh- The cAh branch suffers from large standard model and little higgs background. 
On the other hand, the eZn branch has a very characteristic final state of two isolated ener- 
getic electrons and two higgs. In addition, the eZn branch also has a large branch when the 
lepton yukawa coupling is above 0.5, which is rather difficult to be probed at the LHC [3]. 
The summary of the Feynman diagrams for the signal processes are shown in Figure [5] 



3 Simulation tools 

Signal and background events are generated using the physsim |^ event generator, based 
on the full hclicity amplitudes including gauge boson decays, calculated using HELAS [7] 
and BASES [8], which properly takes into account the angular distributions of the decay 
products. The effects of initial state radiation, beamstrahlung and beam energy spread are 
also included. In this study we will use unpolarized beams. Once the final state particles 
are produced their four momentum are handed down to Pythia 6.4 |10j for parton showering 
and hadronization. For final state tau leptons the decay is carried out by TAUOLA [3] . The 
detector response is simulated using JSFQuickSimulator [IT] ; a fast Monte-Carlo detector 
simulator which implements ILD geometry |l2j and other detector-performance related pa- 
rameters. The assumed detector resolution is specified in Table [21 
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Table 2: ILD detector resolution 
Detector Resolution 


Coverage 


Vctex detector atp = 5.0 {10.0 / pl3sin^^'^ 9) ^m 


cos 6*1 < 0.93 


Central drift chamber crp^/P-^ = 5 x 10"^ PTiGeV/cy^ 


cos6i| < 0.98 


ECAL as/E = 17%/VE 1% 


|cos6i| < 0.99 



HCAL aE/E = 45%/VE © 2% |cos6l| < 0.99 



4 Analysis strategy 

Three physical quantities will be extracted in this study, i.e. model parameters, LHT particle 
masses and LHT couplings. Since the overall analysis procedure is fairly similar among all 
processes, for simplicity we will name the pair-produced LHT parent particle LHp, and 
it's LHT daughter particle LHd. As mentioned before, model parameters can be extracted 
from the measured LHT particle masses |13j . which can be derived by solving kinematics. 
Couplings can be extracted from cross section and angular distribution of LHp. Here we will 
discuss the analysis strategy on how to extract these physical quantities from observables. 

4.1 Mass measurement and parameter extraction 

Due to T-parity conservation, T-parity odd particles are produced in pairs. Once T-Parity 
odd particles are produced(LHp), for the same reason, they decay into a standard model 
particle and a T-parity odd particle(LHd). In general, we can extract both LHT particle 
masses (LHp and LHd) by analyzing the energy distribution of the reconstructed standard 
model particle(SMd). The edges, i.e. the upper and lower kinematic limits of the standard 
model energy {Emax and Emm) are written in terms of LHT particle masses. 

Emax = IhHpEcsM ~^ "lHp^lHpPsm \^) 

Emin = ^LHpE ^^j — PlHp^lHpPsm V"^/ 

where ^lhp^Plhp) is the 7(/3) factor of the LHp in the lab frame, while E*^^{p*^^) is the 
energy (momentum) of the standard model particle in the rest frame of LHp. Keep in mind 
that the standard model particle energy satisfies, ^* = (M^ +AP +M'^ )/(2M, „) 

^ ^-^ ' SMd ^ LHp SMd LHd^' ^ LHp/ 

With this method we can extract new particle masses without any assumption of the model. 
Because the LHT model is described with only two model parameters it gives a characteristic 
mass hierarchy. Performing this model independent mass measurement will provide strong 
evidence that the discovered new particles are indeed LHT particles. 

4.2 Cross section and angular distribution of LHp production 

By assuming the vertex structure, i.e. the spin and the ratio of a right-hand, left-hand 
couplings, we can extract the couplings concerning heavy gauge bosons and heavy leptons 
through cross section measurement. There are a total of eight heavy gauge boson/heavy 
lepton vertacies concerning the five processes we are treating. Therefore in addition to cross 
section measurement we will use the angular distribution of LHp (differential cross section) 
for three processes in order to extract all couplings. 
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The procedure on how to extract the coupUngs are summarized in Figure [3] The red 
circles represent the vertex that will be measured at that process mode and the blue circles 
represent the vertex that are measured in previous process modes. Through differential 
cross section measurement, we can extract information about the couplings contributing to 
the production of LHp, i.e. the coupling sign, whether the contribution is destructive or 
constructive, and the ratio between red circled couplings and the blue circled couplings. 

First, by measuring the cross section of ZhZh process mode, we can extract the coupling 

Oi geenZa- 

Once we have measured this coupling, then we can use it as an input variable to extract 
the geeuAn coupling in the ZhAh process mode. In the cheh process mode we use the 
previously measured couplings geenZH ^-^'^ geenAn ^-s input variables to extract the couplings 
of g Zen en ^-iid 5i/ejfWjf Through differential cross section measurement we can obtain 
the couplings' sign and the coupling ratio and between gzenen ^^^^ other couplings that 
contribute to eue-H production. Now that we have obtained all the couplings used in enSH 
production, we can extract the couplings in the branch, i.e. geHi^Wn through cross section 
measurement 

Coupling extraction in WrWh and vhi^h process modes arc fairly the same. Except, 
for WhWh, cross section measurement is used for measuring the coupling contributing to 
the WhWh production. For i^hi^h, the cross section measurement is used for the extraction 
of the total width of i/h- 

With this procedure, we are able to extract all couplings contributing to the five processes 
that we are treating. 

For simplicity, we will use the following x^ to evaluate the measurement accuracy of the 
couplings. 

Where i? is the observable, C is the coupling, meas is the subscript for a measured value, 
th is the subscript for the true value, i is the superscript that indicates the type of observable, 
e.g. cross section of ZhZh etc. and A is used to represent the standard deviation of the 
measured observable. 

Once we have extracted model parameters from masses and verified the consistency 
among all processes, we can say that it is extremely LHT like and make an estimation on 
the coupling values and search for the x^ minimum point around this estimation point. The 
correlation between two couplings will not be discussed in this study. It will be a subject 
once we have verified the consistency of couplings with model parameters independently. 
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Process mode 



Table 3: coupling measurement strategy 

Branch derived vertex 



observable 



6 ^ Zh 

* ^^^V>'N^>^ 



e 1 Zh 



cross section 






e C2 ^H 

• ^K/sysA- 



cross section 




e 




Wh 






cross section, 
differential cross section 



e ^-v^ w„ 



WmX^N^N^ 



6^"C 






cross section, 
differential cross section 




Wh V 



tot 




tot 



cross section, 
differential cross section 
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5 Results from simulation study 

In this section, we present the resuhs of our simulation study. Simulation has been preformed 
at ^/s=500 GeV for AhZh production and y/s~l TeV for other processes, each process is 
studied with an integrated luminocity of 500 fb~^. 

5.1 Signal selection 

The signal process we are treating include multiple jets in the final state. 



VIS 



The clusters in the calorimeters arc combined to form a jet if the two clusters satisfy yij < 
Vcut- Vij is defined as where 9ij is the angle between momenta of two clusters, £'i(j) are their 
energies, and E^s is the total visible energy. For AhZh process, all events are forced to 
have two jets by adjusting ycut- For other processes, all events arc forced to have four jets, 
which originated from a higgs or a W boson. In order to reconstruct the parent particle we 
have used a x'^ function defined as 

^2 = r^MsMi -*^ MsMf/al^,, + r'MsM2 -'^ MsMf/al,^,, (5) 

where SM is the parent particle, i.e. higgs or W boson, '''^°M5jvfi(2) is the invariant mass of 
the first (second) 2-jet system paired as a parent candidate, *^Msm is the true value of the 
parent particle mass, and omsm is the mass resolution for the parent particle. 

5.2 Selection criteria 

Here we will discuss the selection criteria that is performed in each process mode. The 
selection criteria for AhZh and WhWh process is the same as previous studies [13] . 

1. ZhZh process 

The pair production of Zh is generated with a cross section of 98.0 fb at a center of 
mass energy of 1 TeV. Since Zh decays into Ah and a higgs boson, the signature is two 
higgs bosons in the final state. At a mass with 134 GeV, the higgs boson mostly decays 
to b-jets and W bosons. Therefore in this analysis, we have used a 4-jet final state, which 
originates from the two higgs bosons. We have selected events with x^ below 60 to select 
well reconstructed events. We also require signal events to have no isolated leptons, have 
more than 1 b-tag jet in the final state and have acoplanarity above 20 °. The selection 
criteria is summarized in Tabled 

The requirement of b-tagging is the existence of 2 tracks with 3 a displacement from 
the interaction point, where a is the impact parameter resolution. Wc select events which 
have more than one b-tagged jet. After b-tagging, the dominant background becomes ti. 
We investigated acoplanarity (tt— </>), where (p is the angle between two reconstructed higgs 
candidates in the plane perpendicular to the beam axis. Since the signal has large missing 
momentum due to Ah, there is a large acoplanarity, compared to tt. The isolated- lepton 
rejection is preformed inorder to exclude events where the higgs decays into a W boson, 
followed by W's leptonic decay. These events have large missing energy, which causes the 
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kincniatical lower limit (Emin) in the energy distribution to shift to lower energy region. 
Causing degradation in the measurement accuracy of LHT particles masses. 

In this study, we identified isolated leptons events as leptons which satisfy the following 
equation. 



TrackEnergy < 5 + 25/7 x ConeEnergy 



(6) 



Where TrackEnergy is the energy of the isolated lepton candidate and ConeEnergy is the 
energy around {cos9 > 0.90) the isolated lepton candidate. An isolated lepton track has 
almost no energy surrounding it, as oppose to tracks from a jet where they are accompanied 
by other jet tracks. Equation ([6|) is optimized distinguish between leptons originating from 
H — >■ WW*, where there is at least one W decays leptonicly, and leptons originating from 
H ^bb. 

Table 4: event selection of analyzed processes ZhZh 

tJ of events 
Process cross section (fb) jJ of events after all cuts 



ZrZh 



■ AhAhHH 
AnAHqqqq 



98.0 



49760 



18989 



wwz 


63.86 


31933 


227 


vvWW 


14.67 


7336 


86 


WW 


3069 


1947408 


994 


tt 


192.9 


96472 


1749 


zz 


202.2 


101094 


114 


ZH 


17.98 


8989 


47 


WhWh 


108.6 


54343 


87 



2. ene-H process 

The signal we have used for eu pair production, is where both eu decay to eZn, which 
is further followed by Zh — >■ hAn- The cross section of this process is 4.56 fb at a center of 
mass energy of 1 TeV. Since there are two electrons and two higgs bosons in the final state, 
we have chosen 2 electrons and 4 jets as the signal. We furthermore require the two isolated 
electrons to have opposite charge, the mass region of the two reconstructed particles to be 
between 104 to 164 GcV, and have a missing transverse momentum above 50 GeV. 

After selecting reconstructed mass to be between 104 to 164 GcV, eeWW becomes the 
dominant background. eeWW has small missing transverse momentum, as opposed to the 
signal, where there is a large missing transverse momentum due to Ah- Events before and 
after selection criteria are shown in Table [5] [^ The criteria for selecting a isolated electron 
is to select an electron track with an cone energy below 15 GeV, where the angle of the cone 
satisfies cos 9 > 0.90. The criteria is optimized to maximize selection efficiency of electron 
coming from e//, i.e. 84% and minimize the probability of miss identifying electrons coming 
from b-jet, i.e. 1.2%. 



^The value in 
opposite eharge. 



jj events" colomn showrs the number of events after selecting two isolated electrons writh 
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Table 5: event selection of analyzed processes 

ji of events 
Process cross section (fb) jj of events after all cuts 



effCff -> eeZuZn 








-^ eeAnAnqqqq 


3.91 


1295 


604 


tt 


200.7 


12092 


33 


eeWW 


1020 


10205 





thth 


3.32 


93 


12 


ttZ, ttH, evWZ, eeZZ, WWZ, ZZZ 


114.5 


2461 


13 



3. vrvh process 

The signal we have chosen for vh pair production, is a process where both vh decay 
to eWn, which is further followed by Wh —^ WAh- The cross section of this process is 
55.7 fb at a center of mass energy of 1 TeV. Since there are two electrons and two W bosons 
in the final state, we have chosen 2 electrons (one electron and one positron) and 4 jets as 
the signal. We selected events where the reconstructed W boson mass is between 60 and 
100 GeV. 

Table 6: event selection of analyzed processes 

jj of events 
Process cross section (fb) [I of events after all cuts 



-^ eeAnAHqqqq 



26.0 



13000 



9132 



Vth'^tht^hRh, 3.3 + 3.9 3600 185 

tt,ttZ,UH,evWZ,eeWW,ZZZ,WWZ 2482 1241000 118 



5.3 Mass determination 

The inasses of new particles can be obtained by recognizing upper and lower kinematic 
limits of the SMd energy distribution, from which the masses of the new particles can be 



derived. Figure [2] shows the SMd energy spectrum for each process. Figure 2(a) shows the 
higgs energy distribution coming from the ZhZr process on top of SM background. Figure 



|2(b)|and figure 2(c) shows the electron energy distribution in the ch^^h and vrvh process. 



Here the background is subtracted assuming the background distribution can be understood 

completely. 

From WhWh process mode we extract Wh,Ah masses, from the Z^Zh mode Zh,Ah 

masses, from the enen — >■ gZhcZh modes, ch, Zh masses and from vhvh modes vh^ Wh 

masses are extracted. 



Now when extracting masses from the edge of the energy distribution, we solve kine- 
matics by solving a quadratic equation. This quadratic equation gives two solutions. One 
corresponds to a true mass solution, the other is a false solution. In chgh vhvh process 
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(a) ZhZh process; 
Higgs energy 



(b) eu^H process; 
electron energy 



(c) vh^h process; 
electron energy 



Figure 2: energy distribution 



the false mass solution is small enough that it can be excluded from previous experiments. 
However the false solution in WhWh and ZhZh cannot be excluded by this method. Fig- 
ure 12] shows the contour plot of Ah and Zu mass, which can be obtained by higgs energy 
distribution in the ZhZh process mode. The two neighboring islands correspond to the 
two solutions. In order to select the true solution and determine Wh, Zh and Ah masses, 
we performed a simultaneous fit using the fact that the true solution of the Ah mass is 
common among both process modes, WhWh and ZhZh- Figure S] show the contour plot 
of Ah and Zh mass after performing a simultaneous fit. A single solution can be seen. The 
measurement accuracy of all five LHT particles masses obtained by analyzing the energy 
distribution in four process modes and performing a simulation fit is shown in Table [T] 
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Figure 3: contour plot of Zh and Ah mass 
extracted from ZhZh process 



Figure 4: contour plot of Zh and Ah mass extracted 
from a simultaneous fit of ZhZh, WhWh process 
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Tabic 7: mass measurement accuracy 
particle true value measurement accuracy 



Ah 


81.9 GeV 


1.3% 


Wh 


369 GeV 


0.20% 


Zh 


368 GeV 


0.56% 


en 


410 GeV 


0.46% 


VH 


400 GeV 


0.10% 



5.4 Parameter extraction 



parameter True value measurement accuracy 



/ 

K 



580 GeV 
0.5 



0.16 % 
0.01 % 



Figure 5: measurement accuracy 
of model parameters 
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Figure 6: contour plot corresponding to 1- 
and 2-cr deviations from the best fit point in 
the parameter k and /-plane 

Model parameter / is extracted from the measured Ah, Zh and Wh particle masses. 
Model parameter n is extracted from the measured en masses, which is described by / and 
K. The contour plot of k, and / extracted in the enen process mode is shown in Figure [6l 
The measurement accuracy of extracted model parameters are shown in Table [5j 

5.5 Cross section measurement 

The result of cross section measurement for all five process modes and the measurement 
accuracy for extracted couplings are shown in Table [S] 



11 



LCWSll 



AhZh 


geenAn 


ZrZh 


QeenZH 


eHSH 


QeHiy.WH 


^hvh 
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Table 8: cross section measurement 

process mode extracted vertex cross section meas. accuracy coupling meas. accuracy 

7?70% 3^90% 

0.859% 0.219% 

2.72% 1.40% 

0.949% 0.648% 

0.401% 0.174% 

5.6 Plan. Angular distribution measurement for LHp pair production 

We will measure production angle distribution of LHp in order to measure other couplings. 
To derive the production angles, a quadratic equation is solved using the masses of LHp and 
the momenta of the SM particles, i.e. daughter of LHp with the assumption of a back-to- 
back ejection of the LHp pair. The equation gives either two solutions which contain one 
correct production angle or no solutions when the discriminant of the equation is negative. 
The unphysical negative discriminant comes from misreconstructing SM particle momenta 
or imperfect back-to-back condition of the two LHp particle mainly due to initial state 
radiation. This two dimensional production angle distribution, composed of two solutions 
are compared to those of templates, which are the ideal distribution of production angles 
with different coupling values. We will use x^ (Equation [3|) to derive the x^ minimum point 
and evaluate accuracy. 

6 summary 

The Littlest Higgs Model with T-parity(LHT) is one of the attractive candidates of physics 
beyond the Standard Model since it solves both the little hierarchy and dark matter problems 
simultaneously. One of the important predictions of the model is the existence of new heavy 
gauge bosons and heavy leptons, where they acquire mass terms through the breaking of 
global symmetry necessarily imposed on the model. The determination of the masses are, 
hence, quite important to test the model. Once LHT masses are extracted and it is confirmed 
that they are consistent with the extracted model parameters; global vacuum expectation 
value / and lepton yukawa coupling k, which are the parameters that hold information of 
mass generation mechanism of the LHT model, the next important step will be to measure 
the couplings that appear in the process modes. 

The couplings are described with the two model parameters. Therefore we can verify 
whether the measured couplings are consistent with the couplings estimated from model 
parameters. 

We have performed Monte Carlo simulations in order to evaluate measurement accuracy 
of the heavy gauge bosons and heavy leptons' masses, model parameters and couplings 
at the ILC. The measurement accuracy of LHT particle masses are shown in Table [T] The 
measurement accuracy of LHT model parameters are shown in Table[Sl and the measurement 
accuracy of LHT model coupling through cross section measurement arc shown in Table |8l 

Our plan is to evaluate the measurement accuracy of the other couplings through mea- 
surements of angular distribution. 
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